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ABSTRACT 


LHS 1140 b is the second-closest temperate transiting planet to the Earth with an equilibrium tem- 
perature low enough to support surface liquid water. At 1.730+0.025R&, LHS 1140 b falls within 
the radius valley separating Ho5-rich mini-Neptunes from rocky super-Earths. Recent mass and radius 
revisions indicate a bulk density significantly lower than expected for an Earth-like rocky interior, 
suggesting that LHS 1140 b could either be a mini-Neptune with a small envelope of hydrogen (~0.1% 
by mass) or a water world (9-19% water by mass). Atmospheric characterization through transmission 
spectroscopy can readily discern between these two scenarios. Here, we present two JWST/NIRISS 
transit observations of LHS 1140 b, one of which captures a serendipitous transit of LHS 1140c. The 
combined transmission spectrum of LHS 1140 b shows a telltale spectral signature of unocculted facu- 
lae (5.80), covering ~20% of the visible stellar surface. Besides faculae, our spectral retrieval analysis 
reveals tentative evidence of residual spectral features, best-fit by Rayleigh scattering from an N»- 
dominated atmosphere (2.30), irrespective of the consideration of atmospheric hazes. We also show 
through Global Climate Models (GCM) that H»-rich atmospheres of various compositions (100x, 300x, 
1000xsolar metallicity) are ruled out to >100. The GCM calculations predict that water clouds form 
below the transit photosphere, limiting their impact on transmission data. Our observations suggest 
that LHS 1140 b is either airless or, more likely, surrounded by an atmosphere with a high mean molec- 
ular weight. Our tentative evidence of an Ns-rich atmosphere provides strong motivation for future 
transmission spectroscopy observations of LHS 1140 b. 


Corresponding author: Charles Cadieux 


charles.cadieux. 1G umontreal.ca 


1. INTRODUCTION 


Whether temperate rocky planets orbiting low-mass 
stars have an atmosphere is arguably one of the most 
important scientific questions of the James Webb Space 
Telescope (JWST) mission (Gardner et al. 2023), one 
that can only be answered if a significant fraction of 
its lifetime is dedicated to this endeavour (Cowan et al. 
2015; de Wit et al. 2023; Doyon 2024). Answering this 
question is an essential first step in assessing the hab- 
itability of nearby worlds and identifying the best tar- 
gets for biosignature searches. The first 18 months of 
JWST observations have highlighted both its power and 
versatility for studying the atmospheres of small exo- 
planets. While the first atmospheric reconnaissance of 
TRAPPIST-1 b (Greene et al. 2023) and TRAPPIST-1c 
(Zieba et al. 2023) through eclipse photometry suggests 
that those two planets may be airless worlds, the pres- 
ence of atmospheres on these planets has not been defini- 
tively ruled out (Zieba et al. 2023; Ih et al. 2023; Lin- 
cowski et al. 2023; Turbet et al. 2023; Ducrot et al. 2023). 
The first JWST transmission spectra of the TRAPPIST- 
1 planets, in general, show strong signs of stellar ac- 
tivity in the form of continuum variability, unocculted 
spots/faculae (Lim et al. 2023), and flares (Howard et al. 
2023). Unocculted spots/faculae, responsible for the 
transit light source (TLS) effect (Rackham et al. 2018), 
pose a significant challenge as they can mimic genuine 
atmospheric signals, greatly complicating the interpre- 
tation of transmission spectra. The TLS effect is a re- 
current theme of numerous Cycle 1 and Cycle 2 trans- 
mission spectroscopy programs (e.g., Lim et al. 2023; 
Moran et al. 2023; May et al. 2023; Fournier- Tondreau 
et al. 2024); a fundamental problem that has yet to be 
solved. 

LHS 1140 (Dittmann et al. 2017; Ment et al. 2019; 
Lillo-Box et al. 2020; Cadieux et al. 2024) is a key- 
stone system for habitability studies. This M4.5 dwarf 
(Dittmann et al. 2017) hosts two small transiting plan- 
ets: an outer planet LHS 1140b (Rp = 1.73 Ra), a rare 
object in the middle of the radius valley (Fulton et al. 
2017, 2018; Cloutier & Menou 2020) on a 24.7-day tem- 
perate orbit, and an inner planet LHS 1140c, a warm 
super-Earth (Rp = 1.27 Rẹ) orbiting every 3.78 days. 
The two planets receive 0.43 and 5.3 times the irradia- 
tion of Earth, respectively, with LHS 1140 b comfortably 
situated within the Water Condensation Zone (Turbet 
et al. 2023). Cadieux et al. (2024) recently presented 
a joint study of almost all transit and radial velocity 
observations to date, obtaining a precision of 396 on 
the masses and 296 on the radii of LHS 1140 b and c. 
Among the temperate rocky exoplanets, only those in 
the TRAPPIST-1 system are currently characterized to 


such precision. Their analysis of the composition of 
LHS 1140 b shows that the planet could be enveloped in 
hydrogen (70.196 H/He by mass), akin to a smaller ver- 
sion of the temperate mini-Neptune K2-18 b, which has a 
transmission spectrum characterized by prominent CH4 
and COs absorption bands in the near-infrared (Mad- 
husudhan et al. 2023). Alternatively, LHS 1140 b could 
be a water world with a surface layer of condensed water 
representing 9-19% of the total mass. A 3D Global Cli- 
mate Model (GCM) of the water world scenario predicts 
liquid water at the substellar point for a large range of 
atmospheric compositions (Cadieux et al. 2024). The 
GCM also predicts a distinctive transmission spectrum 
featuring relatively small (~15 ppm) CO» features at 
72.8 and ~4.3 um. 

Transmission spectroscopy with HST/WFC3 shows a 
tentative signal of H50/CH, near 1.4 um in a low mean 
molecular weight atmosphere (Edwards et al. 2021; Bi- 
agini et al. 2024), but a combination of unocculted spots 
and faculae offers an alternative explanation — a hy- 
pothesis also supported by ground-based measurements 
(Diamond-Lowe et al. 2020). These observations imply 
a large fractional coverage (78096) of stellar active re- 
gions which is somewhat surprising given that LHS 1140 
is a relatively old (75 Gyrs; Cadieux et al. 2024), slowly 
rotating star (Pot = 13145 days; Dittmann et al. 2017), 
with low levels of flaring activity (Medina et al. 2022). 
Recent observations of two transits of LHS 1140 b with 
JWST/NIRSpec (G235H/G395H, 1.7-5.2 um) ruled out 
an Ha-rich atmosphere and favored a high mean molec- 
ular weight atmosphere (Damiano et al. 2024). Further 
'TLS characterization with JWST is required, especially 
in the context of seeking the detection of small atmo- 
spheric signals that could be severely affected/biased by 
stellar activity. 

This letter presents the results of a transmission 
spectroscopy program with JWST NIRISS/SOSS (0.6- 
2.8 um) aimed at both characterizing the stellar activity 
of LHS 1140 and providing a strong discriminating test 
of the mini-Neptune/water-world scenarios of Cadieux 
et al. (2024). We show that TLS is clearly detected 
in LHS 1140 and that a cloud-free mini-Neptune atmo- 
sphere for LHS 1140 b is strongly excluded. 

The letter is structured as follows. The observations 
and data reduction steps are presented in Section 2 fol- 
lowed by the transit modeling in Section 3 and spectral 
extraction analysis in Section 4. The retrieved atmo- 
spheric properties of LHS 1140b and TLS constraints 
are presented and discussed in Section 5 followed by 
a discussion of proposed follow-up observations in Sec- 
tion 6. We then conclude in Section 7. 


2. OBSERVATIONS 


'Two transits of LHS 1140 b were obtained with JWST 
using the NIRISS instrument (Doyon et al. 2023) on 
UT2023-12-01 and UT2023-12-26 as part of program 
DD6543 (PI: Cadieux & Doyon). The spectroscopic 
time series were carried using the SOSS mode (Albert 
et al. 2023), offering a spectral resolution R ~ 650, with 
the SUBSTRIP256 subarray configuration to access the 
first two diffraction orders of NIRISS (0.6-2.8 um). Each 
visit comprised 949 integrations, covering 2.58 hr before 
ingress, the 2.15-hour transit and 1.11hr after egress. 
The second visit fully captured a transit of LHS 1140c 
lasting 1.13 hr and starting approximately 34 min after 
first contact of LHS 1140 b. Each integration was com- 
posed of Neroup = 3 filling approximately 50% of the full 
well to minimize charge migration and non-linearity ef- 
fects (Albert et al. 2023). We made use of two indepen- 
dent data reduction pipelines, namely SOSSISSE (Lim 
et al. 2023) and supreme-SPOON (Feinstein et al. 2023; 
Radica et al. 2023, 2024; Benneke et al. 2024), both 
yielding consistent transmission spectra with a mean 
(median) absolute deviation of 0.720 (0.600). More de- 
tails of both pipelines and their comparison are provided 
in Appendix A. The SOSSISSE data products are here- 
after used in the analysis. 

During the first visit, the star LHS 1140 was erro- 
neously positioned outside the target acquisition field 
of view, causing a displacement of the orders on the 
detector with offsets Az — —157px (dispersion) and 
Ay — —12px (cross-dispersion) relative to their nom- 
inal positions. This offset was calculated by compar- 
ing the pixel-to-wavelength calibrations derived from 
the two visits. À modified wavelength solution for the 
shifted spectral trace was determined through a cross- 
correlation between the stellar trace and a PHOENIX 
stellar atmosphere model (Husser et al. 2013), follow- 
ing the procedure described by Feinstein et al. (2023); 
Radica et al. (2023). The target acquisition failure is 
attributable to an incorrect epoch of proper motion en- 
tered during scheduling (LHS 1140 is a high proper mo- 
tion target), which means this sequence was executed 
with blind guiding from the observatory. Fortunately, 
this data set can be fully recovered with a small wave- 
length shift in order 1 — 0.70-2.67 um compared to 
0.85-2.83 um during visit 2 — and negligible change in 
order 2. 

The first visit is also affected by a background star! 
that overlaps the spectral trace of LHS 1140 for wave- 


l Gaia DR3 2371032989200665984, G= 16.667, d= 274.08 pc, 


(Gaia Collaboration et al. 2023) 
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lengths below 0.90 um in order 1, and between 0.60- 
0.75 um in order 2, contaminating the flux at the 1-2% 
level. We estimate a spectral type M2.4+0.4V for the 
contaminant star using the G — Grp calibration of Ki- 
man et al. (2019). The Ter estimate of 3516K from 
Gaia DR3 (Gaia Collaboration et al. 2023) is also con- 
sistent with an M2V spectral type (Pecaut & Mamajek 
2013). We constructed a model of the trace of the con- 
taminant star from a MOV template (WASP-80; GTO 
1201, PI: Lafreniere) observed with NIRISS. This tem- 
plate was aligned by maximizing the correlation coef- 
ficient between the derivatives in spatial direction x of 
the template and observations in regions (pixels) where 
only the contaminant is seen. Once the optimal shift 
was determined, we computed the amplitude of the con- 
tamination pixel by pixel, and use the median of these 
amplitudes to scale the template to match the contami- 
nant trace. We then subtracted this contaminant model 
from the images leaving flux residuals of the order of 
0.196 on the affected pixels of LHS 1140. Such residuals 
could decrease (dilute) the transit depth of LHS 1140 b 
by 10 ppm at most. This possible systematics is further 
mitigated in the next Sections by analysing both vis- 
its together, as visit 2 was not affected by a dispersed 
contaminant. 


3. WHITE LIGHT CURVE ANALYSIS 


One output of SOSSISSE is an amplitude time se- 
ries a(t) that corresponds to broadband flux measure- 
ments (details in Appendix A), hereafter referred to as 
the ‘white’ light curve (WLC). The WLC is modeled 
with two components, one for transits and the other 
for treating systematic signals, both described below. 
The fit was performed using the juliet framework (Es- 
pinoza et al. 2019) that generates transit models with 
batman (Kreidberg 2015) and implements Gaussian Pro- 
cess (GP) regression with celerite (Foreman-Mackey 
et al. 2017). The posterior distribution is explored with 
nested sampling and 500 live points using the dynesty 
package (Speagle 2020) available in juliet. 

The transit component assumes circular orbits for 
both planets (Gomes & Ferraz-Mello 2020). The orbital 
parameters of planet k (k: ‘b’, ‘c’) are the period Py 
fixed to the value from Cadieux et al. (2024), the time 
of inferior conjunction to,,, and the scaled semi-major 
axis a,/R, that we derive using a common stellar den- 
sity p,. For p,, we use the same Gaussian prior as in 
Cadieux et al. (2024) constructed from stellar mass and 
radius estimates. For all other parameters, we adopt 
wide uniform priors. The transits of LHS 1140b and 
c are parameterized with the planet-to-star radius ra- 
tio Ry ,/R, and the impact parameter by. We fit for 
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quadratic limb-darkening parameters qı and q2 defined 
in Kipping (2013) that ensure physical solutions for val- 
ues between 0 and 1. 

For the systematic component, we use GP regression 
to model two nuisance signals also reported in other 
NIRISS data sets (e.g., Coulombe et al. 2023; Radica 
et al. 2024; Benneke et al. 2024): a beat pattern in- 
troduced by the observatory thermal control (McElwain 
et al. 2023) and correlated structures on the timescale 
of a few minutes suggestive of stellar granulation (e.g., 
Kallinger et al. 2014; Grunblatt et al. 2017; Pereira et al. 
2019). We use a combination of two simple harmonic 
oscillator (SHO) kernels in celerite to jointly model 
these systematic signals. We fix the period of the first 
SHO term to 204s, corresponding to the dominant peri- 
odicity in the out-of-transit fluxes based on its power 
(highest peak) in the generalized Lomb-Scargle peri- 
odogram (Zechmeister & Kürster 2009). We fit for an 
amplitude opeat and a quality factor Qpeat for this os- 
cillation. As in Radica et al. (2024), we use a critically 
damped SHO term (Q fixed to V/0.5) for the granulation- 
like signal and fit for an amplitude ogra and a timescale 
Tgra Of this stochastic variation. An extra jitter term 
jitter is added in quadrature to the diagonal of the co- 
variance matrix to account for excess noise. We also fit 
for a baseline flux and a temporal slope for each visit. 
We measure slopes of —0.26 and —0.36 ppm per min 
during visits 1 and 2, respectively, a flux variation ap- 
proximately 20-30 times weaker than observed during 
the transits of TRAPPIST-1 b with NIRISS (Lim et al. 
2023). This is consistent with the fact that LHS 1140 is 
a slow-rotating star with a flux level remarkably stable 
over a timescale of a few hours, as also confirmed by 
TESS (Ricker et al. 2015). 

'The WLC of the two visits were jointly fit with a com- 
mon px, q1, and q2 to allow more precise constraints. 
The adopted priors and resulting posteriors (16^, 50th, 
and 84t}? percentiles) of this fit are presented in Table 1. 
The best-fit transit and systematic components of the 
WLC are shown in the top panels of Figure 1. Fit- 
ting visit 1 and visit 2 independently yielded consistent 
parameters within the lo uncertainty, but we note a 
systematic increase of ~40 ppm in the transit depth of 
LHS 1140 b during visit 2 due to the double transit. The 
joint fit alleviates the covariance between Rp», Rp,c, and 
b. in visit 2 by using the information of a single transit of 
planet b from visit 1. This same argument is later used 
to extract combined transmission spectra of LHS 1140 b 
and c by jointly fitting the two visits in each spectral 
channel. 


Table 1. Transit and systematic parameters inferred 
from a joint fit of the white light curve of visits 1 and 2 


Parameter Prior Posterior 

Stellar parameters 

ps (g: cm?) N (26.0, 2.6?) 26.0 5$ 
Orbital parameters 

P, (days) (fixed)* 24.73723 

P, (days) (fixed)* 3.777940 

to,» (BJD-2460300) U (4.65, 4.75) 4.69824(3) 

to,c (BJD-2460300) U (4.65, 4.75) 4.70775(4) 

ay / Ry (derived)! 5.2") 

ac/ R. (derived)! DEI 


Transit parameters 


U(0,50) 7.439 + 0.015 
U(0,50) 5.312 + 0.028 


Ry,» / R+ (96) 
Rp,c/Rx (96) 


bp U (0, 0.9) git. 
be U (0,0.9) Do 
qı u (0,1) 0.16 + 0.02 
q2 U (0, 1) 0.22 + 0.05 
GP parameters 
beat (ppm) LU (0.1, 1000) 69 +6 
beat LU (0.1, 1000) 113727 
gra (ppm) LU(0.1, 1000) 5345 
Tera (min) LU(1, 100) 5.2 + 1.0 
Tjitter (ppm) LU(0.1, 1000) 95 +3 


Note—*Period from Cadieux et al. (2024). ‘Scaled 
semi-major axis (a/R,) derived from p, and P 


4. SPECTROPHOTOMETRIC ANALYSIS 


We binned the extracted spectroscopic time series 
from SOSSISSE to R ~ 100 in both order 1 and order 
2. We excluded wavelengths below (above) 0.86 um in 
spectral order 1 (2), as these wavelengths are already 
covered by the other order at a higher throughput. Be- 
cause of the wavelength coverage differing between the 
two data sets (see Sect. 2), this results in 145 spec- 
trophotometric light curves for the first visit and 149 
for the second visit. We also constructed a common 
wavelength grid at the native resolution corresponding 
to the average of the nearest pixels (in wavelength) be- 
tween the two visits. This average grid deviates at most 
by 0.23nm from the wavelength calibrations from indi- 
vidual visits. We combined the spectral time series of 
visit 1 and visit 2 using the common wavelength grid, 
then binned to a R œ% 100 to end up with 142 joint 
spectrophotometric light curves. We then removed the 
systematic component determined from the WLC from 
all spectral light curves. At R ~ 100, the median pho- 
tometric precision of the light curves is approximately 
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Figure 1. Spectrophotometric transit fits of LHS 1140b for visit 1 (UT2023-12-01, left) and visit 2 (UT2023-12-26, right) and 
LHS 1140c for visit 2 only. The broadband (‘white’) light curves are shown in top panels with the best-fit full model (transit 
+ systematics) shown in black and systematics component only in pink (scaled 10x for clarity). For both visits, examples of 10 
systematics-corrected spectroscopic bins are depicted with colored points in bottom panels. 


800 ppm which prevents us from characterizing the beat 
pattern and granulation signals, both with amplitudes 
below 100 ppm (Table 1). 

'To derive the transmission spectra of LHS 1140 b and 
c, we fitted a transit model to each spectrophotomet- 
ric light curve. Again, the fits were carried out with 
juliet using nested sampling and 500 live points for 
each spectral bin. We fixed Pk, tor, ar/Rx, bj to the 
best-fit values from the WLC analysis (Table 1) and fit 
for Ry. / R,. The limb-darkening coefficients were calcu- 
lated for all spectral channels with ExoTic-LD (Grant & 
Wakeford 2022) using a custom PHOENIX model input 
(Tee = 3100 K, log g = 5, [Fe/H] = 0.0) for LHS 1140. 
For the spectral fits, the corresponding qı and q» limb- 
darkening parameters are kept fixed. We allow for a 
baseline flux parameter and a Cjitter for each visit and all 
spectroscopic bins. Examples of 10 spectroscopic chan- 
nels jointly fitting the transits of LHS 1140 b and c are 
shown in Figure 1. 

'The final combined transit spectrum of LHS 1140b 
is presented in Figure 2 along with a binned version 


at R ~ 20 to ease comparison with HST/WCF3 data. 
The spectrum is relatively flat above 1 um, but features 
a clear 200-ppm decrease towards shorter wavelengths. 
We show in the next section that this is a signature of the 
TLS effect from unocculted faculae. The HST/WFC3 
data are fairly consistent with NIRISS albeit for one or 
two data points with larger transit depth blueward of 
1.4um which may have largely contributed to the ten- 
tative water/methane detection mentioned by Edwards 
et al. (2021) and Biagini et al. (2024). The transmission 
spectra of LHS 1140 b and c are available in Appendix B. 


5. THE ATMOSPHERE OF LHS 1140B 


We now use the combined transmission spectrum of 
LHS 1140 b to jointly infer properties of the planetary 
atmosphere and stellar contamination from unocculted 
active regions resulting in stellar contamination (the 
transit light source, or ‘TLS’ effect). Our analysis con- 
siders the evidence for stellar contamination (Sect. 5.1), 
evidence for a planetary atmosphere through retrievals 
(Sect. 5.2), and finally a comparison with GCM pre- 
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dictions of mini-Neptunes (Sect. 5.3). We also provide 
a preliminary analysis of the transmission spectrum of 
LHS 1140c in Appendix C showing that its spectrum 
is featureless and inconsistent with a cloud-free Hs-rich 
atmosphere. 


5.1. Transit Light Source Effect from LHS 1140 


We first consider whether LHS 1140 b's transmission 
spectrum can be explained only by unocculted stel- 
lar active regions. We modeled the TLS effect us- 
ing the transmission spectra stellar contamination mod- 
ule of the POSEIDON retrieval package (MacDonald & 
Madhusudhan 2017; MacDonald 2023). "We generate 
model spectra by multiplying a bare-rock transmission 
term, (Ry/R,)?, by the wavelength-dependent stellar 
contamination factor from two distinct stellar hetero- 
geneities (Equation 3, Fournier-Tondreau et al. 2024). 
The active region spectra are calculated by interpolat- 
ing PHOENIX stellar models using the PyMSG package 
(Townsend & Lopez 2023). Our TLS model fits for the 
covering fractions and temperatures of the facula, spot, 
and quiet photosphere, yielding six free parameters: Rp, 
frac, Fspot; Trac; Tspot, and Tphot- We fix the surface grav- 
ity of all active regions to logg = 5.04 (Cadieux et al. 
2024). The priors are listed in Table D1 and mainly 
follow those in Fournier-Tondreau et al. (2024) for their 
two-heterogeneity model. Our POSEIDON model spectra 
are generated at R = 20000 from 0.60 to 2.9 um, con- 
volved with a Gaussian kernel to the native resolution 
of NIRISS/SOSS (R = 650), then multiplied by the in- 
strument sensitivity function and binned down to the 
wavelength spacing of the data (R z 100). The pa- 
rameter space is explored via the nested sampling pack- 
age PyMultiNest (Feroz & Hobson 2008; Buchner et al. 
2014) with 2000 live points. The posterior distributions 
for this TLS-only model are shown in orange in the bot- 
tom row of Figure 2 (also in Appendix D), with the 
median and 1c confidence intervals listed in Table D1. 

LHS 1140 b's transmission spectrum shows definitive 
evidence of unocculted stellar faculae. We detect facu- 
lae at 5.80 confidence (Aln Z = 14.91, where A In Z is 
the increase in log Bayesian evidence for stellar contam- 
ination over a flat line). Our 6-parameter TLS model 
(x? = 1.73, with 136 degrees of freedom), shown in Fig- 
ure 2, provides a significantly better fit compared to a 
flat line (x2 = 1.96, with 141 degrees of freedom). The 
reduced x? of 1.73 may be indicative of underestimated 
error bars, or perhaps imperfect TLS modeling due to 
inherent uncertainties in stellar atmosphere models (Lim 
et al. 2023; Jahandar et al. 2024). We find unocculted 
faculae 72183 K hotter than the photosphere covering 
20+17% of the visible stellar disk. There is no evidence 


for spots. A wide range of covering fractions are consis- 
tent with the data due to the degeneracy between ffac 
and AT = Trac — Tpnot (i.e., smaller /larger faculae need 
to be hotter/cooler to produce a similar spectral con- 
tamination). This level of stellar activity is consistent, 
within a factor of order unity, with the observed pho- 
tometric peak-to-peak variability of ~1% of LHS 1140 
(Dittmann et al. 2017) inferred from a 2-year monitoring 
campaign with MEarth (Irwin et al. 2009). 

We note that fitting a single stellar contamination 
model to LHS 1140b’s combined spectrum amounts to 
assuming a similar surface distribution of unocculted 
heterogeneities during the times of the two visits sep- 
arated by one orbital period of 24.7 days, equivalent 
to ~20% of one stellar rotation. To validate this as- 
sumption, we also fit the same TLS-only model to the 
transmission spectrum of each individual visit, obtain- 
ing fic = 0.181018 and Tk, = 315170 K for visit 1 
and fme = 0.06'018 and Tm. = 3420722: K for visit 
2. We also observe consistent unocculted faculae prop- 
erties from LHS 1140c's transmission spectrum with 
fu. 012 075 and Trac 820327 K during visit 2 
(Appendix C). This confirms that our two visits have 
consistent faculae properties, and hence it is valid to an- 
alyze the combined spectrum with a single TLS model. 
The consistent faculae parameters from our LHS 1140 b 
and LHS 1140c analyses confirm that these unocculted 
faculae lie outside the stellar surface sampled by the 
two planetary transit chords. We additionally note that 
Damiano et al. (2024) reported a spot crossing event 
during a different set of two transits of LHS 1140 b ob- 
served with NIRSpec in July 2023. Their observations, 
conducted 124 days prior to our NIRISS data set (about 
one full rotation of LHS 1140), disfavor unocculted het- 
erogeneities (their Table 4). However, the longer wave- 
lengths covered by NIRSpec G235H and G395H are less 
sensitive to the TLS effect. Our analysis therefore crit- 
ically underscore the need for NIRISS observations to 
characterize the nature and spectral impact of unoc- 
culted stellar active regions on transmission spectra. 


5.2. Atmospheric Retrieval 


We next establish whether there is any evidence of 
an atmosphere on LHS 1140 b. Given the presence of 
unocculted faculae established above, we performed a 
joint stellar contamination and planetary atmosphere re- 
trieval analysis using POSEIDON. We consider a range of 
atmospheric retrieval models, including single-gas pure 
atmospheres, multiple gases, and atmospheres with or 
without aerosols. Our models consider the following 
gases (opacity sources in parentheses): CO» (Tashkun 
& Perevalov 2011), H20 (Polyansky et al. 2018), CH4 
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Figure 2. Combined transmission spectrum of LHS 1140b with NIRISS/SOSS from two transits and POSEIDON retrieval 
analysis results. The binned spectrum at R ~ 20 (white points) and at a higher R ~ 100 (gray points) significantly improve 
upon the HST/WFC3 data from Edwards et al. (2021) (black diamonds). The median models with 1-20 confidence envelopes 
are shown for a stellar contamination model (orange) and a joint stellar contamination and planetary atmosphere model (blue), 
with their respective posterior distributions shown in the bottom panels (same colors). Confidence intervals for each parameter 
are given in Table D1. The transmission spectrum of LHS 1140b is mainly shaped by stellar contamination from unocculted 
faculae. A planetary atmosphere dominated by Na or CO» is compatible with the data, with the former pure composition 
being the maximum a posteriori model. A clear N2-rich atmosphere combined with unocculted faculae is preferred at 2.30 over 


faculae alone (see Sect. 5.2 and Fig. D1). 


(Yurchenko et al. 2017), Hə (Hohm 1994 for Rayleigh 
scattering; Karman et al. 2019 for collision-induced ab- 
sorption), and Na (Sneep & Ubachs 2005 for Rayleigh 
scattering; Karman et al. 2019 for collision-induced ab- 
sorption). For multi-gas models, the volume mixing 
ratios follow a permutation-invariant centered-log-ratio 
(CLR) prior (Benneke & Seager 2012), which allows for 
any molecule to be equally likely a priori to be the dom- 
inant gas. We assume isothermal atmospheres and fit 
for a reference radius, Rp yer, at a reference pressure of 
l bar. We include a two-parameter power-law prescrip- 
tion for hazes (MacDonald & Madhusudhan 2017) in all 
multi-gas atmospheric models. For the cloudy model, 
we assume an optically-thick gray opacity/surface with 
all layers deeper than Pj set to infinite opacity. We 
allow Mp to vary as a free parameter following a Gaus- 


sian prior (5.60 +0.19Mẹ; Cadieux et al. 2024). For the 
stellar contamination, we use the same TLS parameter- 
ization as outlined above (i.e., 5 parameters to fit for 
both faculae and spots). We used 2000 PyMultiNest 
live points for all retrievals. A summary of the free pa- 
rameters and their priors is provided in Table D1. 

Our retrieval analysis identifies tentative evidence for 
a cloud-free Na-dominated atmosphere on LHS 1140 b 
(2.30). We first considered models with a single-gas 
atmosphere with negligible cloud opacity in the upper 
atmosphere (all including stellar contamination), find- 
ing that the fit did not improve for CO2, H20, CH3, 
or Hə pure atmospheres. For example, our 100% COs 
model has a lower Bayesian evidence compared to stel- 
lar contamination only (In Z — 1141.02 vs. 1142.15) — 
reflecting the ‘Occam penalty’ from adding one or more 
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parameters that do not improve the fit (see e.g., Trotta 
2008) — since there are no apparent CO» features in our 
NIRISS spectrum. Pure CO2 atmospheres must have 
T « 233K (20 upper limit) to be compatible with our 
non-detection of CO» absorption (i.e., high-temperature 
pure CO» atmospheres are ruled out). However, our 
100% Nə model has the highest Bayesian evidence of 
all our models (In Z = 1144.65), with a Aln Z = 2.50 
(2.80) for No + stellar contamination over stellar con- 
tamination alone. The evidence for Nə arises from a 
combination of a residual short-wavelength slope com- 
patible with Na Rayleigh scattering and a weak fea- 
ture near 2.2 um attributable to Na-Na collision-induced 
absorption (CIA) — see Appendix D and Figure D1. 
Our multiple gas retrieval reaffirms that Na is the fa- 
vored background gas, even when hazes are included 
in the model, with a retrieved abundance of 100*996 
(see Fig. 2, blue posteriors), indicating that our data 
rule out the tail to lower Na abundances present in 
the CLR prior. We rule out cloud-free H2O-rich, CH4- 
rich, and Hs-rich atmospheres, with 20 upper limits of 
log H20 < —2.94, log CH4 < —2.78, and log Hə < —0.95 
(« 1196), respectively. While our multiple gas retrieval 
does allow for cold (T ~ 100K) COs-rich atmospheres 
(corresponding to a relatively featureless atmospheric 
contribution to the spectrum), this solution is not sta- 
tistically favored, as discussed above. To provide a more 
conservative estimate for the Na detection significance, 
we calculated the Bayesian evidence for a multiple gas 
model including all the other gases and hazes but with- 
out Na (ensuring the possibility of both CO2-dominated 
atmospheres and hazes are propagated into the Bayesian 
evidence), finding a conservative A ln Z = 1.50 (2.30) in 
favor of an N5-dominated atmosphere. Our tantalizing 
evidence for a N5-dominated atmosphere on a habitable 
zone super-Earth raises a resounding call for additional 
observations of LHS 1140 b to confirm this result. 

We additionally considered the potential for aerosols 
in LHS 1140 b's atmosphere. First, we reiterate that 
our multi-gas retrieval discussed above included a para- 
metric prescription for atmospheric hazes in the form of 
a power law (MacDonald & Madhusudhan 2017). The 
haze parameters are unconstrained in all our retrievals 
(see posteriors in Appendix D). The reason for this is 
that hazes still exist within an atmosphere, providing 
an additional opacity enhancing the dominant Rayleigh 
slope, but the retrieval must first ‘select’ a background 
gas before the impact of a haze can be established. If 
Na is the background gas, Na Rayleigh scattering alone 
provides a good fit, and additional haze opacity is not 
needed (hazes are redundant parameters). If CO» is the 
dominant gas, then hazes have no spectral effect due to 


the high molecular weight and low temperature shrink- 
ing the scale height and hence flattening the spectrum. 
Therefore, allowing for hazes does not affect the infer- 
ence of a clear No-dominated atmosphere as the favored 
model. When we include a cloud/surface in our multi- 
gas retrieval we obtain a bimodal solution: (i) a clear 
^ 10096 Na atmosphere, as above, or (ii) an optically- 
thick cloud deck at low Puf ^ 107° bar with no con- 
straints on the atmospheric composition (i.e., such high- 
altitude clouds would render any transmission spectrum 
featureless, regardless of the atmospheric composition). 
'The Bayesian evidence decreases when adding a cloud 
deck (see Table D1), indicating that this is a redun- 
dant parameter not required to explain LHS 1140 b's 
transmission spectrum. We also stress that such ex- 
ceedingly high-altitude clouds are probably unrealistic. 
According to 3D GCMs of LHS 1140b for both the 
water world (Cadieux et al. 2024) and hydrogen-rich 
scenarios (see below), a cloud deck (e.g., H20, CO2 
clouds) should form much deeper in the atmosphere 
(P = 1bar). Therefore, while our present observations 
do not strongly favor a clear No-dominated atmosphere 
over a high-altitude cloud deck, the latter scenario is less 
plausible from a GCM standpoint. 

We examine the evidence for Na in the atmosphere of 
LHS 1140b and present our full atmospheric retrieval 
results and posterior distributions in Appendix D. 


5.3. 9D GCM Forward Modeling 


Here we compare the observed spectrum with real- 
istic predictions from 3-D GCMs of a mini-Neptune 
with a 80bar Hə-rich atmosphere of various composi- 
tions: 100x, 300x, and 1000xsolar metallicity. For this, 
we used the Generic PCM (Planetary Climate Model), 
historically known as the Generic LMD GOM, as in 
Cadieux et al. (2024). The model is capable of simu- 
lating various types of exoplanets, including terrestrial 
planets in the TRAPPIST-1 system (Turbet et al. 2018; 
Fauchez et al. 2019) and sub-Neptunes such as GJ 1214 b 
(Charnay et al. 2015) or K2-18b (Charnay et al. 2021). 
'The GCM simulations performed here closely follow the 
methodology described in Charnay et al. (2021), except 
that star and planet properties have been adapted to 
LHS 1140b. The GCM experiments are described in 
more details in Appendix E. 

Following the methodology of Fauchez et al. (2019), 
the outputs of the GCMs are used as inputs to the Plane- 
tary Spectrum Generator (Villanueva et al. 2018) to gen- 
erate synthetic transmission spectra using the GlobES 
module (Villanueva et al. 2022) that consider the 3D 
nature of the atmosphere and cloud in a self-consistent 
manner. The simulations include atmospheric refrac- 
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Figure 3. Stellar contamination-corrected combined spectrum of LHS 1140 b using the best-fit solution of Figure 2. The black 
points are binned to R ~ 20 while the gray points are at a higher R ~ 100. Dashed lines depict GCM-based forward models 
for mini-Neptune atmosphere (100x, 300x, 1000xsolar metallicity) and solid lines the water world GCMs of Cadieux et al. 
(2024) for an Earth-like atmosphere (1 bar N2, 400 ppm CO2) and a pure CO» atmosphere at 5 bar. An H>-rich atmosphere for 
LHS 1140 b is formally rejected (7106) by the NIRISS/SOSS data. A flat solution (e.g., water world or airless planet) provides 
a significantly better fit but the reduced x? of 1.75 may indicate an underestimation of error bars, an imperfect TLS correction, 
or the presence of residual spectral features (e.g., slope). An apparent linear trend with a slope of approximately 40 ppm across 
the NIRISS domain (dotted green line) improve the fit (x? = 1.72), but is not predicted by our GCMs. 


tion, collision-induced absorption (CIA) and multiple 
scattering. Irrespective of atmospheric composition, all 
forward models have a cold trap between ~0.1 and 1 bar 
(see Fig. El), favoring condensation of H20 into clouds. 
For these realistic cloudy mini-Neptune models, CH4 
and CO» (to a lesser extent) are the only observable 
transmission features in the near-infrared. Figure 3 
compares the TLS-corrected transmission spectrum of 
LHS 1140b with these updated mini-Neptune mod- 
els along with water world GCMs taken from Cadieux 
et al. (2024) for an Earth-like atmosphere (1bar Na, 
400 ppm CO2) and a pure CO» atmosphere at 5 bar. All 
hydrogen-rich GCMs are formally rejected by the data 
(fitting for an offset) with a confidence level exceeding 
100, with a Aln Z = —326.90, —172.22, —84.98 for the 
100x, 300x, and 1000xsolar metallicity models, respec- 
tively, compared to a flat spectrum. The water world 
GCMs or a flat line (airless planet) provide a much bet- 
ter fit to the data. However, we note an apparent slope 
of ~40 ppm in the spectrum that is not predicted by the 
GCMs explored in this study. The possible origin of this 
residual slope is discussed in Appendix D. 


6. DISCUSSION 
6.1. A Secondary Atmosphere on LHS 11400? 


While an exceptionally cloudy mini-Neptune is com- 
patible with our NIRISS spectrum, such a scenario 
would be in stark contrast with the relatively clear 
atmospheres predicted by self-consistent GCMs of 
LHS 1140b. Recent observations of the temperate 


mini-Neptunes K2-18b (Madhusudhan et al. 2023) and 
TOI-270d (Benneke et al. 2024) have detected strong 
spectral features associated with Ho-rich atmospheres, 
with relatively little cloud opacity. Instead, our tenta- 
tive inference of an Nj-dominated atmosphere suggests 
LHS 1140b is a water world with a secondary atmo- 
sphere (Forget & Leconte 2014; Kite & Ford 2018; Kite 
& Barnett 2020; Marounina & Rogers 2020). 

One line of evidence supporting atmospheric reten- 
tion on LHS 1140b is that the planet is likely resilient 
to atmospheric mass loss of light elements (Ha, He) from 
processes such as XUV-driven photoevaporation (Owen 
& Wu 2017) and core-powered mass loss (Ginzburg et al. 
2018). According to Cadieux et al. 2024 (Figure 4, 
therein), an initial envelope mass fraction as low as 
~0.1% H/He could mostly persist after 10 Gyr thanks to 
the relatively low instellation of LHS 1140 b combined 
with its relatively high surface gravity. Based on these 
results, it is probable that an atmosphere composed of 
any elements heavier than hydrogen or helium would 
also be retained by the planet over similar timescales. 

Another way to present this argument is illustrated 
in Figure 4, showing the instellation vs. the escape 
velocity of LHS 1140b (S = 0.43 + 0.0888, ve = 
20.13 +0.37 km/s; Cadieux et al. 2024) along with other 
exoplanets smaller than 1.8 Rẹ using the NASA Exo- 
planet Archive (Akeson et al. 2013). This radius cut 
is chosen to focus on terrestrials, super-Earths, and po- 
tentially water worlds, excluding exoplanets predomi- 
nantly composed of gas. The empirical “cosmic shore- 
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line" (S o vå.) of Zahnle & Catling (2017) based on 
Solar System objects appears as an orange band in the 
figure, a region that separates planets with an atmo- 
sphere from those without. Reconnaissance atmospheric 
characterization in transmission and/or in emission with 
JWST of a few key targets will be required to establish 
this shoreline across spectral types. As shown in this 
diagram, LHS 1140 b, TOI-1452 b (Cadieux et al. 2022), 
TRAPPIST-1f and TRAPPIST-1 g (Agol et al. 2021) 
appear as the most favorable small exoplanets to host 
an (secondary) atmosphere. 
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Figure 4. Escape velocity-instellation diagram of small ex- 
oplanets (Rp < 1.8R@). The empirical “cosmic shoreline” 
(S x vi.) of Zahnle & Catling (2017) based on Solar Sys- 
tem bodies is depicted in orange with an arbitrary envelope 
to underscore uncertainty with spectral types. The size of 
data points scale with their Transmission and Emission Spec- 
troscopy Metrics (Kempton et al. 2018) to highlight favor- 
able targets for JWST atmospheric characterization. Targets 
with S < 108, or with TSM/ESM above 90" percentile 
are annotated. The TRAPPIST-1 system is abbreviated as 
“T-1”. The optimistic habitable zone (Early Mars/Recent 
Venus) defined by Kopparapu et al. (2013) is shown in green 
(approximated for a range of spectral types). With a rela- 
tively high escape velocity and low instellation, LHS 1140 b 
likely retained an atmosphere. 


6.2. Future observations 


Considering the significant implications associated 
with the potential detection of a secondary atmosphere 
on LHS 1140 b, seeking unambiguous evidence warrants 
a combined approach of eclipse and transmission spec- 


troscopy. Here, we discuss the required observations to 
ascertain with high certainty the presence of a secondary 
atmosphere on LHS 1140 b (or lack thereof). It should 
be noted that due to the system's limited visibility with 
JWST, only 4 or 5 transits/eclipses of LHS 1140 b are 
observable in a given year. 


6.2.1. Eclipse Photometry 


Eclipse photometry centered on the 15-um CO» ab- 
sorption feature is a powerful tool for constraining the 
presence of an atmosphere on small rocky planets, e.g., 
TRAPPIST-1b (Greene et al. 2023) and TRAPPIST- 
lc (Zieba et al. 2023). The eclipse depth is given by 
Fy/F, where F, and F, are respectively the flux from 
the dayside of the planet and the star during the eclipse. 
Assuming a BT-COND stellar atmosphere model (Al- 
lard et al. 2012) with an effective temperature of 3100 
K anchored to the 2MASS J magnitude of LHS 1140, 
we find F} = 10.1 mJy for the F1500W MIRI filter. 
The observations of TRAPPIST-1b at 15 um suggest 
that the stellar model prediction for F, should be ac- 
curate to „10% (Greene et al. 2023). A blackbody 
is assumed for Fp with the equilibrium temperature 
Toq = Tog \/Rs/a [f(1— Ag)] ^, where Tag is the ef- 
fective temperature of the host star, R, its radius, a 
the planet semi-major axis, Ap the planet Bond albedo 
and f the re-radiation factor that can take two extreme 
values: f — 1/4 corresponding to a homogeneous dis- 
tribution of the energy across the planet and f — 2/3 
corresponding to no heat redistribution to the night side 
(Esteves et al. 2013). 

Figure 5 shows the predicted emission spectrum of 
LHS 1140 b for two airless models with different albedos 
and the two water world GCM cases of Cadieux et al. 
(2024). Three visits with MIRI would be sufficient to 
detect/exclude at the 30 level the most optimistic sce- 
nario (Ap = 0, f = 2/3) corresponding to a dark, airless 
planet which is the coreless case inferred from the inter- 
nal structure model of Cadieux et al. (2024). An airless 
Europa-like surface (Ag = 0.6, Buratti & Veverka 1983) 
would require ^19 visits to detect the secondary eclipse, 
i.e., observing every eclipse for nearly 5 years. Detecting 
and discriminating the water world cases (an Earth-like 
or COs-rich atmosphere) is out of reach through eclipse 
photometry at 15 um only, and would probably require 
the combination of many filters (e.g., 15 + 21 jum) over 
many eclipses. 


6.2.2. "Transmission spectroscopy 


Cadieux et al. (2024) estimated that 12 transits (six 
times each of NIRSpec G235 and G395) would be re- 
quired to detect the CO» features longward of 2.8 um, 
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Figure 5. Predicted emission spectrum of LHS 1140b for 
various scenarios including the Earth-like (1bar N2, 400 
ppm CO2) and pure CO; (5 bar) atmosphere cases from the 
GCMs of Cadieux et al. (2024). Extreme Bond albedos cases 
(0 or Europa-like) with no atmospheric heat redistribution 
(f — 2/3) are shown with dashed lines. The grey enveloppe 
represents the 1096 uncertainty of LHS 1140's stellar flux in- 
ferred from a BT-COND model. The red cross, centered on 
the deepest airless case (Ap — 0, f — 2/3), gives the predicted 
uncertainty at 15 wm for three visits needed to reach a 3-0 de- 
tection assuming photon-limited performance. Similarly, the 
green cross gives the number of visits needed (seven) for de- 
tecting the deepest eclipse at 21 jum. Each visit is assumed to 
last 7.5 hours (1hr + 3 times the transit duration), a simi- 
lar experimental design adopted by Greene et al. (2023) and 
Zieba et al. (2023) for the MIRI observations of TRAPPIST- 
1b and c, respectively. 


an exposure time in line with the nine additional G395 
visits proposed by Damiano et al. 2024 to detect the 
CO» feature at 4.3 um. While G395 may be the most 
attractive mode for maximizing the signal-to-noise of 
the 4.3 um feature, it comes with the significant scien- 
tific compromise that the atmospheric water abundance 
in the presence of TLS, as observed in LHS 1140 b, is 
likely to be poorly constrained. The same ambiguity is- 
sue between TLS and water detection associated with 
the G395 mode alone has been noted by other NIRSpec 
programs on warm super-Parths (e.g., Moran et al. 2023; 
May et al. 2023). Moreover, the G395 mode is insensi- 
tive to short wavelength spectral slopes associated with 
hazes and Rayleigh scattering. Our observations pro- 
vide tantalizing evidence (2.30) for a Rayleigh scatter- 
ing slope and a CIA feature around 2.2 jum that could 
be explained by a Ns-rich atmosphere. Such a signal 
could be confirmed with approximately 4 additional vis- 
its that should also include NIRSpec observations since 
the No-N2 CIA shows stronger opacity at 4.3 um which 
is very close to the CO» feature. Clearly, the optimal 
experimental design for exploring the full diversity of 
secondary atmospheres (N2- or COs-rich) along with 
proper TLS characterization calls for the widest pos- 
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sible wavelength coverage as offered by the combina- 
tion of NIRISS/SOSS and NIRSpec/G395. Observing 
within the same season (July or December) through two 
consecutive visits between NIRISS and NIRSpec mini- 
mizes the change in stellar rotation phase which is only 
20% given the relative orbital period of LHS 1140 b 
(24.7 days) and the stellar rotation period of 131 days. 
'The two consecutive transits of LHS 11140 b observed 
with NIRISS that we presented in this letter is the proof 
of concept of this experimental design, showing that the 
TLS can vary weakly within one observing season, im- 
plying that the G395 and NIRISS data acquired within 
one season could be analyzed jointly under the reason- 
able assumption that both visits share a similar stellar 
heterogeneity configuration. 


7. SUMMARY & CONCLUSION 


We have presented the 0.65-2.7 um transmission spec- 
trum of the temperate planet LHS 1140 b obtained from 
two visits with JWST/NIRISS, the second of which coin- 
cidentally involved a transit of LHS 1140c. All spectra 
exhibit a low level of stellar contamination caused by 
unocculted faculae with covering fractions at the level 
of ~20%. GCMs of mini-Neptunes with various atmo- 
spheric compositions (100x, 300x, 1000xsolar metal- 
licity) are all excluded with a significance greater than 
100. A spectral retrieval analysis also excludes a clear 
H>-rich atmosphere, with the most likely atmospheric 
scenario being that of an Nə- or COs- dominated en- 
velope. The former Ns-rich atmospheric composition is 
favored by the data at 2.30 from a tentative detection of 
Na Rayleigh scattering combined with weak Na-Na CIA 
absorption near 2.2 um. While an Ha-rich atmosphere 
with a relatively high cloud deck is consistent with the 
NIRISS observations, such a muted spectrum is in con- 
tradiction with cloud formation simulated in GCMs and 
the relatively clear atmospheres detected on the temper- 
ate mini-Neptunes K2-18b (Madhusudhan et al. 2023) 
and TOI-270d (Benneke et al. 2024). 

The NIRISS observations are more likely consistent 
with the water world scenarios presented by Cadieux 
et al. (2024), with or without a secondary atmosphere, 
a conclusion also supported by recent transmission spec- 
troscopy data obtained with JWST/NIRSpec (Damiano 
et al. 2024). Of all nearby small temperate planets, 
LHS 1140b is the most likely to have retained a sec- 
ondary atmosphere based on its low instellation and 
high surface gravity. The next obvious step to better 
constrain LHS 1140 b's atmospheric composition is to 
perform a joint analysis of both NIRISS and NIRSpec 
datasets, ideally with a common data reduction method- 


ology. 
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Only three visits of 15 um eclipse photometry with 
MIRI is required to yield a 3-0 detection of the sec- 
ondary eclipse associated with a dark, airless planet, 
but nearly five years (4—5 eclipses per year) is required if 
its albedo is closer to an icy surface (Ag ~ 0.5) as sug- 
gested by the relatively high water mass fraction (15%) 
of LHS 1140b. Transmission spectroscopy is the most 
efficient method to detect the potential secondary at- 
mosphere of LHS 1140 b through consecutive observa- 
tions between NIRISS/SOSS and NIRSpec/G395 dur- 
ing the same season to allow for a proper characteriza- 
tion of stellar contamination at all epochs and explore 
Nə- and COs-dominated atmospheres. Given the lim- 
ited visibility of LHS 1140 b, several years worth of ob- 
servations may be required to detect its potential sec- 
ondary atmosphere, a program that should be initiated 
as soon as possible given the limited lifetime of JWST. 
LHS 1140b is arguably the best temperate transiting 
planet for which liquid surface water may be indirectly 
inferred through the detection of sufficient level of at- 
mospheric CO». 
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APPENDIX 


A. DATA REDUCTION PIPELINES 


The NIRISS data were reduced with two independent pipelines (SOSSISSE and supreme-SPOON) briefly described 
here. 

The SOSSISSE package starts with the rateints.fits data product issued by the Mikulski Archive for Space Telescopes 
(MAST) and performs custom background subtraction, flat-field correction, and 1/f noise correction. A high signal- 
to-noise, normalized reference trace M; j is then constructed from stacking the out-of-transit corrected 2D images, 
with i,j the pixels in the x (dispersion) and y (cross-dispersion) directions. Here, Mi; is limited to 30 pixels in 
the y-axis. In the SOSSISSE framework, flux is expressed as a linear combination of M; j and its spatial derivatives 
tracing instrumental morphological changes of the trace (Equation Al, Lim et al. 2023), e.g., x- and y-displacement, 
rotation, PSF width changes. We utilized two outputs from SOSSISSE: the amplitude term a(t), scaling M; j at each 
integration equivalent to the (normalized) stellar flux integrated over all pixels (‘white’ light curve), and the extracted 
2D spectroscopic time series for the spectrophotometric fits. 

The supreme-SPOON pipeline works with the uncal.fits files available on MAST. We followed Stages 1, 2, and 3 
detailed in Radica et al. (2023) and extracted the flux using a simple box aperture with a width of 30 pixels. The 
aperture box size was selected to minimize the variance in the out-of-transit white light curve. We expect a negligible 
level of order contamination (order 2 onto order 1) at wavelengths above 1.6 um for this target (Darveau-Bernier et al. 
2022; Radica et al. 2022). A 1/f noise correction is applied following the method described in Radica et al. (2023); 
Benneke et al. (2024). We, moreover, correct for the presence of the dispersed contaminant star in visit 1 using the 
method of Radica et al. (2023). This correction is done after flux extraction, unlike the SOSSISSE treatment of the 
contaminant which is applied directly on the detector images (Sect. 2). 

Figure Al shows that the combined transmission spectra at R ~ 100 derived from both pipelines are in good 
agreement with mean (median) absolute deviations of 0.72 ø (0.60 0) for LHS 1140 b and 0.580 (0.520) for LHS 1140c. 
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Figure Al. Combined (two visits) transmission spectrum of LHS 1140b on the left and (one visit) transmission spec- 
trum of LHS 1140c on the right extracted with the SOSSISSE and supreme-SPOON pipelines here showed for the first order 
of NIRISS/SOSS. The transit light curve fitting is identical between pipelines; any differences arise during the reduction stages. 
The panels below show the difference in transit depths (SOSSISSE — supreme-SPOON) normalized by the maximum errors between 
the two pipelines at each spectral bin. A running mean (window = 5) is drawn with a blue curve showing no clear trends in 
the residuals. The pipelines are in good agreement with an average per-point difference of 0.72 and 0.580 for the spectra of 
LHS 1140b and c, respectively. 
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Table B1. NIRISS SOSS transmission spectra of LHS 1140b and c 


Central wavelength (um) Bandwidth (um) Depth b (ppm) Error b (ppm) Depth c (ppm) Error c (ppm) Order 


0.65338 0.00327 5210.71 236.64 3331.44 362.37 2 
0.65995 0.00330 5592.47 253.80 2247.57 397.11 2 
0.66658 0.00333 5094.17 227.03 3275.11 391.74 2 
0.67328 0.00337 5167.24 199.42 2676.46 374.10 2 
0.68005 0.00340 5296.80 175.59 2650.32 358.80 2 
0.68688 0.00343 5112.36 165.52 2645.80 321.08 2 
0.69379 0.00347 5322.05 141.06 2833.85 259.10 2 


NoTE— Table B1 is available in its entirety in machine-readable format. 


B. TRANSMISSION SPECTRA 


The final NIRISS SOSS transmission spectra (R ~ 100) of LHS 1140 b and c from the combined spectral extraction 
are given in Table B1. 


C. PRELIMINARY ANALYSIS OF THE TRANSMISSION DATA OF LHS 1140C 


'The first transmission spectrum of the warm super-Earth LHS 1140c is presented in Figure C1. We present here 
a first exploratory analysis of LHS 1140c’s NIRISS/SOSS spectrum by repeating the retrieval exercise outlined in 
Section 5.2 for the Flat line, TLS-only, and TLS + Multi-gas (Haze) models. All priors remain unchanged except for 
Rp,ref and Mp changed to U (1,1.5) and N (1.91, 0.062) based on Cadieux et al. (2024), and T' now allowed to go to 
higher temperatures (100—600 K). 

A flat solution is favored for LHS 1140c (2.10) due to the lower signal-to-noise ratio spectrum derived from a 
single transit entirely confined into a transit of LHS 1140 bb (double transit, Fig. 1). The Flat line model yields a 
In Z = 1090.49, while the TLS-only and TLS + Multi-gas (Haze) retrieval runs end up with marginally lower Bayesian 
evidences of In Z — 1089.34 and 1088.83, respectively. 
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Figure C1. Transmission spectrum of LHS 1140 c from NIRISS/SOSS at a R ~ 20 (white points) and at a higher R ~ 100 (gray 
points). The best-fit models from POSEIDON with 1-20 confidence envelopes are shown for a stellar contamination model (orange) 
and a joint stellar contamination and planetary atmosphere model (blue). The TESS and Spitzer photometric measurements 
(Cadieux et al. 2024) are shown in red. We have included the Spitzer point in the retrieval analysis, but its exclusion does not 
change the results. The NIRISS/SOSS observations is flat with marginal evidence of unocculted faculae consistent with the 
simultaneous transit of LHS 1140 b (Fig. 2). The spectrum covering 0.65-2.7 um has an average transit depth that falls between 
those measured by TESS and Spitzer, resolving the 4c discrepancy in radius measurements between the two instruments noted 
by Cadieux et al. (2024). 
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The stellar contamination only model yields consistent unocculted heterogeneity properties compared to LHS 1140 b: 
frac = 0.12*00 and Trac = 3293779? K, with fspor consistent with 0. The TLS effect associated with LHS 1140c is 
smaller by definition given the surface ratio of the planets (R),c/Rp,»)” of about one half. 

Our joint stellar contamination and planetary atmosphere retrieval (TLS + Multi-gas model) can inform on possible 
atmospheres on LHS 1140c. The posterior distributions for this run are presented in Appendix D. We do not detect 
any molecular feature or a slope from hazes (haze parameters are unconstrained). The most likely composition favored 
by the data is a pure Ng envelope at any T. Two other solutions exist that are approximately equally likely: a H2O-rich 
atmosphere at T ~ 200 K or a COs-rich atmosphere at T < 500 K. The combination of a high mean molecular weight u 
and/or low T' for such atmospheres is sufficient to flatten the spectrum below the sensitivity of our NIRISS spectrum. 
The data is incompatible with clear CH4-rich and H2-rich atmospheres with 20 upper limits of log CH4 < —2.25 and 
log Hə < —0.94 (< 11%), respectively. 


D. SUPPLEMENTARY MATERIAL: ATMOSPHERIC INFERENCE ANALYSIS 
D.1. Evidence for a N5-Dominated Atmosphere on LHS 114006 


Here, we elucidate why our retrievals favor an Ng-dominated atmosphere. The observed transmission spectrum 
is given by A) obs = €x, star 4d, atm, Where €, star is the ‘stellar contamination factor’ (e.g., Equation 3 in Fournier- 
Tondreau et al. 2024) and A) atm is the ‘regular’ transmission spectrum of the planetary atmosphere. Given that 
unocculted faculae are the dominant wavelength-dependent effect sculpting LHS 1140 b's transmission spectrum, we 
first corrected our NIRISS/SOSS data by dividing out the best-fitting unocculted stellar contamination model (binned 
to the resolution of the data). This allows one to compare model atmosphere-transmission-only spectra, A, atm, to 
an equivalent form of the data. 

Figure D1 compares our stellar contamination-corrected transmission spectrum of LHS 1140 b to our best-fitting 
100% Nə and 100% CO, atmosphere models. The CO» model is essentially flat, consistent with a non-detection 
of COs bands. The Nə model, however, exhibits a clear Rayleigh slope in both the first and second orders of the 
NIRISS/SOSS data (alongside a weaker N3-N» collision-induced absorption feature near 2.2 um). Both models allow 
for the possibility of power-law hazes, but Nə Rayleigh scattering alone provides the necessary slope without any 
additional scattering opacity. We stress that the residual slope attributed here to Nə relies on the assumption that 
the stellar contamination is well-described by our interpolated PHOENIX models for faculae. The best TLS model 
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Figure D1. Tentative evidence of an N?-dominated atmosphere on LHS 1140 b. The stellar contamination corrected LHS 1140 b 
NIRISS/SOSS transmission spectrum (black points for R ~ 20, gray points for R ~ 100) is compared to the best-fitting model 
transmission spectra for a 100% Nz atmosphere (green) and a 100% CO» atmosphere (blue), plotted at R = 20. The N2 model 
has a best-fitting T = 399 K, while the CO2 model has T = 116 K. The Na model provides a better fit (2.30) due to the presence 
of a Rayleigh scattering slope and a weak N5-N» collision-induced absorption feature. 
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leaves a significant residual near 2.3 um which is coincident with the CO bandhead, a well-known temperature- and 
gravity-sensitive spectral feature in M dwarf spectra. Given inherent uncertainties associated with stellar atmosphere 
models (e.g., Lim et al. 2023; Jahandar et al. 2024), one cannot rule out the possibility that the measured slope is an 
artifact of imperfect stellar contamination correction. The evidence for Na is only tentative at this stage (2.30) and 
will need to be confirmed with future observations. 


D.2. Atmospheric and Stellar Retrieval Results 


The priors, posterior parameter constraints, and key statistics from our retrieval analysis of LHS 1140 b are presented 
in Table D1. The full posterior distributions are shown for the TLS-only model in Figure D2, TLS + Nə (Clear) in 
Figure D3, TLS + Multi-gas (Haze) in Figure D4, and TLS + Multi-gas (Haze + Cloud) in Figure D5. Finally, the 
posterior distributions for the TLS + Multi-gas (Haze) retrieval of LHS 1140c are presented in Figure D6. 


Table D1. Retrieval analysis summary for LHS 1140 b's combined NIRISS/SOSS transmission spectrum 


Parameter Prior Flat line 'TLS only TLS + N2 TLS + Multigas TLS + Multi-gas 
(Clear) (Haze) (Haze + Cloud) 


Stellar Contamination 


Tspot (K) U (2300, Ter + 307,4) - 3003468, 2995+69 3017+69 3010+59 
Tphot (K) N (Tof, 07,4) - 3073+35 3074+33 3076+33 3073+34 
Trac (K) U (Teg — 3 OTs, 1-2 Toft) = 3155725 3153793 3146795 3142+66 
Jspot U (0,0.5) = 0.04* 013 om 0.07405 0.06005 
fiac U (0, 0.5) = 0.20+9:17 0.22+9-15 0.25+9-15 es aa 


Atmospheric Properties 


Rp,ret (Re) U (1.5, 2.0) 1.7513+0.0006 1.762+0.002 1.765+0.002 1.764+0.002 1.758+0.005 
M (Me) N (5.60, 0.19?) - - 5.59+0.17 5.48+0.17 5.60+0.16 
TI en) - s 331153 29075, 224 
log nase U (74,8) : g - E 2 
Yhaze U (—20, 2) = = — —8.35+6-52 —8.34*832 
log(P;urr/bar) U (-7,2) u = = _ Sein) 
log Xu, (derived) = = - « —0.95 « —0.56 
log Xx, CLR (—12, 0) - - 0 (fixed) 0:00 -0.69* 0-59 
log Xco, CLR (—12,0) = = — -7.621373 1 garter 
log Xu;o CLR (712,0) = - - < —2.94 -6.144347 
log Xon, CLR (—12,0) = = = < —2.78 -6.171519 


Model Statistics 


Nparams 1 6 8 14 15 
x? (d.o.£.) 277 (141) 235 (136) 228 (134) 228 (128) 228 (127) 
x 1.96 1.73 1.70 1.78 1.80 
In Z 1127.24 1142.15 1144.65 1141.79 1141.60 


NOTE—TLS = ‘Transit Light Source Effect’ (contamination from unocculted stellar active regions, Rackham et al. 2018). *Multi- 
gas’ retrievals include N2, Ha, CO», H20, and CH4. CLR refers to the centered-log ratio prior (e.g., Benneke & Seager 2012). 
Hə also follows a CLR prior, even though it is not a free parameter (it is derived from the summation to unity condition), since 
this is defining feature of the centered-log ratio. The a priori known stellar properties are Tog = 3096 and or, = 48K (Cadieux 
et al. 2024). + CO» has a biomodal distribution, with either (i) a non-detection for a pure N2 atmosphere or (ii) 100% CO» 
to render the spectrum flat (the latter solution is statistically disfavored). Upper bounds correspond to 95% confidence (2c) 
limits. Figure 2 shows results for the ‘TLS only’ (orange histograms) and the ‘TLS + Multi-gas (Haze)’ (blue histograms) models 
summarized in this Table. 


20 


= +0.002 
Rp, ret = L762 00 


EIE 


É oot 


fg, = 0.207037 
T T T 


| 
i} 
l 
1 
| 
1 


fac 


Trot = 3003715, 


T spot 


Ta, = 3155435 


Trac 


T phot = 3073435 


T phot 


SS © 49 
y y ku y 


Rp, ref fspot fiac 


Figure D2. Posterior distribution for the stellar contamination only (‘TLS-only’) POSEIDON retrieval of LHS 1140 b's combined 
NIRISS/SOSS transmission spectrum. 
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Figure D3. Posterior distribution for the joint Na atmosphere and stellar contamination retrieval of LHS 1140 b's combined 
NIRISS/SOSS transmission spectrum. Note that No does not show in the corner plot as a free parameter because the Na 
abundance is fixed to 10096 for this model. 
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Figure D4. Posterior distribution for the joint multi-gas atmosphere with hazes and stellar contamination retrieval of 
LHS 1140 b's combined NIRISS/SOSS transmission spectrum. 
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Figure D5. Posterior distribution for the joint multi-gas atmosphere with hazes, a cloud/surface, and stellar contamination 
retrieval of LHS 1140 b’s combined NIRISS/SOSS transmission spectrum. 
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Figure D6. Posterior distribution for the joint multi-gas atmosphere with hazes and stellar contamination of LHS 1140c’s 
NIRISS/SOSS transmission spectrum. 
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E. GLOBAL CLIMATE MODEL SIMULATIONS AND SYNTHETIC OBSERVABLES 


We performed a series of 3-D GCM (Global Climate Model) simulations designed to explore a large variety of 
plausible atmospheric compositions for LHS 1140 b. These simulations include (i) thick mini-Neptune atmospheres 
(80 bar Ha-rich) with compositions of 100x, 300x and 1000xsolar metallicity, described in the main text, (ii) compact 
secondary Nz and CO»s-rich atmospheres, described in Cadieux et al. (2024), and (iii) sensitivity experiments of H5-rich 
atmospheres on top of a dry surface or a global surface ocean a.k.a. ‘hycean’-type planet (Madhusudhan et al. 2021). 

Figure E1 shows the temperature, water vapor mixing ratio and cloud profiles for the mini-Neptune cases, as well as 
for the Earth-like and CO2-dominated atmospheres for comparison. In all these simulations, cloud decks are located 
below atmospheric pressures of ~0.1 bar, near the top of the tropospheric cold trap. Clouds therefore have a very 
limited impact on the transit spectra (see Fig. 3) for all the simulations explored in this work. Previous intercomparison 
works (Sergeev et al. 2022; Fauchez et al. 2022) have shown that cloud deck altitude can vary slightly from one model 
to another, but here cloud deck altitude would have to rise by at least two orders of magnitude to alter our conclusions. 
The presence of high-altitude hazes — which is not simulated here — could flatten the transit spectrum even in the 
case of an atmosphere rich in H5, but in this scenario we would expect to see a haze slope in the shortwave part of 
spectrum (Sing et al. 2011). 
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Figure E1. Results of Global Climate Model (GCM) simulations of LHS 1140 b assuming mini-Neptune thick Hz atmosphere 
(100x, 300x and 1000xsolar metallicity composition), an Earth-like atmosphere (1 bar N2, 400 ppm CO2), and a thick CO2- 
dominated atmosphere (5 bar of CO2). Panels show vertical profiles of the atmospheric temperatures (left), water vapor mixing 
ratio (middle), and water cloud mixing ratios (right). The thick solid lines indicate the global mean vertical profiles, and the 
dotted lines indicate the terminator vertical profiles (impacting transit spectra). 


Some simulated cases show a slight variability in the amplitude of transit spectra due to cloud variability (in time 
and location), as already shown in Charnay et al. 2021 in the case of K2-18 b. However, even at model timesteps where 
the terminator region is the cloudiest, the amplitude of the transit spectrum features for the mini-Neptune cases is far 
greater than the scatter of NIRISS observations, so these compositions can be formally rejected. 

We ran additional sensitivity simulations in which we changed the boundary conditions for the mini-Neptune simula- 
tions. First, we added a dry surface at 10 bar atmospheric pressure, but in this case the synthetic spectra (not shown) 
are very similar to the mini-Neptune cases (see Fig. 3). Second, we added a surface entirely covered by an ocean 
also at 10 bar atmospheric pressure. This aquaplanet endowed with a thick H2-dominated atmosphere, also known 
as hycean planet, has a quite different behavior. In fact, the simulation enters runaway greenhouse, which forces the 
ocean to evaporate and the surface temperatures to rise. This result is compatible with previous calculations of the 
runaway greenhouse limit for hycean planets (Innes et al. 2023). As the simulation progresses, the amount of water 
vapor in the atmosphere increases at all altitude layers. The altitude of the water cloud deck also increases, similar 
to what has been shown already for No-dominated atmospheres entering the runaway greenhouse (Chaverot et al. 
2023). For all the hycean cases we have simulated, the amplitude of the spectral atmospheric features is far greater 
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than the dispersion of the NIRISS/SOSS spectrum, even accounting for the effect of clouds. However, these runaway 
greenhouse simulations have by definition not reached their final equilibrium state, with clouds potentially forming 
even higher in the atmosphere. Moreover, it is not known whether the ocean at the base of the H5-H2O-dominated 
atmosphere can eventually stabilize at higher temperatures, which would also likely affect the water vapor and cloud 
structures. Last but not least, we have not accounted for convection inhibition, which would likely make this type 
of atmosphere much hotter in the lowest atmospheric layers, and even less likely for an ocean to exist (Innes et al. 
2023; Leconte et al. 2024). Performing fully converged 3-D GCM simulations of hycean planets, accounting for moist 
convection, clouds, and convection inhibition together at the same time, is a challenging yet necessary step (Leconte 
et al. 2024) to simulate accurately the evolution and stability of LHS 1140 b as an hycean planet. 


